Numerical evidence of breaking 
of vortex lines in an ideal fluid 



EVGENIY A. KUZNETSOV 1 - 3 , OLGA M. PODVIGINA 2 ' 3 
& VLADISLAV A. ZHELIGOVSKY 2 ' 3 

1 L.D. Landau Institute for Theoretical Physics, 
2 Kosygin str., 117334 Moscow, Russian Federation 
kuznetso@itp.ac.ru 

o 

International Institute of Earthquake Prediction Theory 
and Mathematical Geophysics, 

79 bldg. 2 Warshavskoe ave., 113556 Moscow, Russian Federation; 
Laboratory of general aerodynamics, Institute of Mechanics, 
Lomonosov Moscow State University, 
1, Michurinsky ave., 119899 Moscow, Russian Federation 

3 Observatoire de la Cote d'Azur, CNRS UMR 6529, 
BP 4229, 06304 Nice Cedex 4, France 

Abstract Emergence of singularity of vorticity at a single point, not related to any 
symmetry of the initial distribution, has been demonstrated numerically for the first 
time. Behavior of the maximum of vorticity near the point of collapse closely follows 
the dependence (to — t) , where to is the time of collapse. This agrees with the 
interpretation of collapse in an ideal incompressible fluid as of the process of vortex 
lines breaking. 

1. Introduction 

The problem of collapse in hydrodynamics, i.e. of a process of sin- 
gularity formation in a finite time, is essential for understanding of the 
physical nature of developed turbulence. Despite a progress in construc- 
tion of statistical theory of Kolmogorov spectra within both diagram and 
functional approaches (see, e.g., Monin & Yaglom 1992; L'vov 1991 and 
references therein), so far the question whether the Kolmogorov spec- 
trum is a solution to the statistical equations of hydrodynamics remains 
open. Another important problem, as yet unsolved, is the one of inter- 
mittency. In statistical sense intermittency can be interpreted as a con- 
sequence of a strongly non-Gaussian distribution of turbulent velocity, 
resulting in deviation of exponents for higher correlation functions from 
their Kolmogorov values (Frisch 1995). Non-Gaussian behavior implies 
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that odd correlation functions do not vanish; this indicates the presence 
of strong correlations between velocity fluctuations, suggesting existence 
of coherent structures in turbulence. Analysis of both numerical and ex- 
perimental data reveals (see Frisch 1995 and references therein) that 
in the regime of fully developed turbulence distribution of vorticity is 
strongly inhomogeneous in space - it is concentrated in relatively small 
regions. What is the reason of this? Can such a high concentration be 
explained by formation of singularity of vorticity in a finite time? How 
can one derive from this hypothesis the Kolmogorov spectrum? This 
question is not rhetoric: it is well known that any singularity results 
in a power-law kind of spectrum in the short-scale region. Thus, the 
problem of collapse is of ultimate importance in hydrodynamics. 

The most popular object in the studies of collapse in hydrodynamics 
is a system of two anti-parallel vortex tubes, inside which vorticity is 
continuously distributed (Kerr 1993), or in a more general setup - flows 
with a higher spatial symmetry (Boratav & Pelz 1994; Pelz 1997). It is 
well known, that two anti-parallel vortex filaments undergo the so-called 
Crow instability (Crow 1970) leading to stretching of vortex filaments in 
the direction normal to the plane of the initial distribution of vortices and 
to reduction of their mutual distance. It was demonstrated in numerical 
experiments (Kerr 1993) that point singularities are formed in cores of 
each vortex tubes at the nonlinear stage of this instability, and \u\ near 
the point of collapse increases like (to — t) -1 , to being the time of collapse 
(see also Grauer, Marliani & Germaschewski 1998). 

2. Basic equations 

In this paper we present results 1 of a numerical experiment, which can 
be interpreted as emergence of singularity of vorticity at a single point in 
a three-dimensional ideal hydrodynamic system, where initial data lacks 
any symmetry. The representation of the Euler equation for vorticity 
u>(r,t) in terms of vortex lines is employed, which was introduced in 
Kuznetsov & Ruban 1998: 

W (r,t) = ( Wo (a)-V a )R(a,t)/J. (1) 

Here the mapping 

r = R(a,t) (2) 

represents transition to a new curvilinear system of coordinates asso- 
ciated with vortex lines, so that b = (u>o(a) • V a )R(a, t) is a tangent 
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vector to a given vortex line, J = det || dR/da. || is the Jacobian of the 
mapping (2). Dynamics of the vector R(a,i) satisfies 

d t K = Uv(K,t), (3) 

where v(R,t) is the flow velocity at a point r = R and II is the transverse 
projection to the vortex line at this point: 

n Q/3 = 5 a(3 - £ = b/|b|. (4) 

Equations (l)-(3) are closed by the relation between vorticity and veloc- 
ity: 

w(r,t) = V x v(r,t); V • v = 0. (5) 

The system of equations (l)-(5) can be regarded as a result of partial 
integration of the Euler equation 

9 ( w = Vx[vxw], V-v = 0. (6) 

A vector field u^o(a) incorporated in (1), V a - u>o(a) = 0, is the Cauchy in- 
variant, manifesting frozenness of vorticity into the fluid. If R(a, 0) = a, 
ujq is the initial distribution of vorticity. 

The Jacobian J can take arbitrary values because the description un- 
der consideration is a mixed, Lagrangian-Eulerian one (Kuznetsov & 
Ruban 1998; Kuznetsov & Ruban 2000). In particular, J can vanish 
at some point, which by virtue of (1) implies a singularity of vortic- 
ity. It was demonstrated by Kuznetsov & Ruban 2000 that collapses of 
this type are possible in the three-dimensional integrable hydrodynamics 
(Kuznetsov & Ruban 1998), where in the Euler equation (6) a modified 
relation between vorticity and velocity (both generalized) is assumed: 

v = Vx(6H/6u>), H = j\u\dY. (7) 

Emergence of singularity of vorticity at a point, where J = 0, means 
that a vortex line touches at this point another vortex line. This is the 
process of breaking of vortex lines. Being analogous to breaking in a gas 
of dust particles (dynamics of a gas with a zero pressure), this process 
is completely determined by the mapping (2). 
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3. Breaking of vortex lines 

Let us assume now that collapse in the Euler hydrodynamics occurs 
due to breaking of vortex lines. Denote by t(a) > a solution to the 
equation J(a, t) = 0, and let to = min a t(a), where the minimum is 
achieved at a = ao- Near the point of the minimum (to, ao) the Jacobian 
can be expanded (cf. Kuznetsov & Ruban 2000): 

J = a (t -t)+ jijAcii Acij + (8) 

where a > 0, 7 is a positive definite matrix and Aa = a— ao. The Taylor 
expansion (8) is obtained under the assumption that J is smooth, which 
is conceivable up to the moment of singularity formation. At t = to the 
numerator in (1), i.e. the vector b, does not vanish: the condition J = 
is satisfied when the three vectors dH/dai (i = 1, 2, 3) lie in a plane, but 
generically none of them equals zero (that were a degeneracy) so that 
near the point of singularity 

b(t ,a ) , . 

<^( r )*)~^ N A A • (9) 

v ' a(to -t) + jijAaiAaj v ; 

Furthermore, J = implies that an eigenvalue of the Jacoby matrix 
(say, Ai) vanishes, and generically the other two eigenvalues (A2 and A3) 
are non-zero. Therefore, there exist one "soft" direction associated with 
Ai, and two "hard" directions associated with A2 and A3. It follows from 
(8), that in the auxiliary a-space the self-similarity Aa ~ (to — t) 1 ^ 2 
is uniform in all directions. However, in the physical space the scales 
are different. Following Kuznetsov & Ruban 2000, we show how an 
anisotropic self-similarity emerges in the flow. The analysis for the Euler 
equation coincides with that for the integrable hydrodynamics (7). 

Decompose the Jacoby matrix J in the bases of eigenvectors of the 
direct (Jip( n ^ >= A n ?/' <n - ) ) and conjugate (ip( n U = Xnip^) spectral prob- 
lems: 

^-^ = E A ^ (n) ^ n) - ( 10 ) 

* n=l 

The two sets of eigenvectors are mutually orthogonal: 

(^W.^M) =S nm . 

In a vicinity of the point of collapse the eigenvectors can be regarded as 
approximately constant. 

Decompose the vectors x and V a in (10) in the respective bases, de- 
noting their components by X n and A n : 

X n = (x-V (n) ), ^- = ^ (n) ' V «)- 
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The vector Aa can be represented in terms of A n as follows: 

Aa Q = ^Vi n) l^ (n) | 2 ^n. 

n 

As a result, (10) can be expressed as 
dX 1 



dA 1 



t + F mn A m A n , (11) 



9X2 -A 9X3 -A H21 

where 

r m „ = 7 Q/3 (A 2 A 3 )-Vi n) ^ m) ^ (n) l 2 l^ (m) l 2 

and r = a(io — t)/(\2\s) is assumed to be small. Consequently, size 
reduction along the directions ip^ and ip^ is the same as in the auxiliary 
a-space, i.e., r 1 / 2 , but in the soft direction, ip^\ the spatial scale is 
~ r 3 / 2 . Therefore, in terms of new self-similar variables £i = Ai/r 3 / 2 , 
(2 = X2/T 1 / 2 and C3 = X3/T 1 / 2 , integration of the system yields for £2 
and C3 a linear dependence on r\ = Aa/r 1 / 2 , and for £i - a cubic one: 

Ci = {l + T ij ri i r) j )ri 1 + -r li r li r l j + -r 11 r)f, i,j = 2,3 (13) 

C2 = A 2 r? 2 , (3 = A3J73. (14) 

Together with (9), relations (13) and (14) implicitly define the depen- 
dence of u on r and t. The presence of two different self-similarities 
shows, that the spatial vorticity distribution becomes strongly flattened 
in the first direction, and a pancake-like structure is formed for t — > to. 
Due to (1) and the degeneracy of the mapping (J = 0), vorticity cj lies 
in the plane of the pancake. Near the singularity behavior of oj is defined 
by the following self-similar asymptotics: 

cj=r- 1 0(Ci,C2,C 3 ). (15) 

In essence, in the above analysis one is concerned with the behavior of the 
mapping near a fold, and thus breaking of vortex lines can be naturally 
explained within the classical catastrophe theory (Arnold 1981; Arnold 
1989). 
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4. Super-weak collapse 

According to the collapse classification of Zakharov & Kuznetsov 1986, 
breaking of vortex lines is not a weak collapse but a super-weak one, 
because already a contribution from the singularity to the enstrophy 
/ = J \u\ 2 dr characterizing the energy dissipation rate due to viscosity is 
small, contribution to the total energy is ~ t ' . However, the 

integral f |Vo;| 2 <ir is divergent as t — > to- Thus, the breaking solution 
v = v(r, t) cannot be continued beyond t = to m the Sobolev space 
H 2 (R 3 ) with the norm || f || ? = (E g <2 / I V<?f ?dv) 1/2 . According to the 
theorem proved by Beale, Kato & Majda 1984, this suffices for 

to 

sup \u\dt = oo (16) 

r 

to hold. The condition (16) is necessary and sufficient for collapse in the 
Euler equation, and it is satisfied for (15). 

Another restriction follows from the theorem by Constantin, Feferman 
& Majda 1996, stating that there is no collapse for any t G [0, to] if 

f °sup|V£| 2 dt < oo, (17) 
Jo 

where the supremum is over a region near the maximum of vorticity 
Occurrence of collapse implies divergence of the integral (17) for 
r — ► 0. Consequently, sup|V£| has to increase at least like r™ 1 / 2 . It 
is evident that, due to solenoidality of u>, either the derivative (£ • V)£ 
in the direction along the vector u> in the pancake-like region should 
have no a singularity at the scales of the order of r 1 / 2 or larger, or the 
singularity should be weaker than r -1 / 2 . However, this does not rule 
out large gradients of £ in a region separated in the soft direction from 
the pancake-like region, for instance, with the behavior d£/dX\ ~ r~ a 
with 1/2 < a < 3/2. This conjecture is plausible, since transition from 
the a-space to the physical one involves a significant contraction in the 
soft direction of the region near the point of breaking: a sphere of radius 
~ r 1 / 2 is mapped into the pancake-like region. Thus, a sphere in the 
r-space of radius ~ r 1 / 2 containing the pancake includes a large preim- 
age of the region outside the sphere in the a-space of radius ~ r 1 / 2 (the 
shape of the preimage is governed by higher order terms in the expansion 
(8) ). Hence in the process of breaking of vortex lines three scales can 
appear: l\ ~ r 3 / 2 , l± ~ r 1 / 2 and an intermediate scale ~ r Q with 
1/2 < a < 3/2 (the presence of which assures that there are no contra- 
dictions with the theorem of Constantin, Feferman & Majda 1996). 
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5. Numerical results 

To verify the hypothesis that formation of singularity in the solutions 
to the Euler equation can be due to vortex line breaking, we performed a 
numerical experiment for the system of equations (1-5). Two features of 
this system are notable. First, in contrast with the original Euler equa- 
tion, possessing an infinite number of integrals of motion - the Cauchy 
invariants, - the system (1-5) is partially integrated and therefore con- 
tains the Cauchy invariants explicitly. Hence, while the invariants are 
guaranteed to be conserved when (1-5) is solved numerically, it is nec- 
essary to test to which extent they are conserved in the course of direct 
numerical integration of the Euler equation (6). Second, in the system 
(1-5) integration in time (in (3) ) is separated from integration over space 
(in (5) ), i.e. from inversion of the operator curl. 

The system (1-5) is considered under the periodicity boundary condi- 
tions and inversion of the operator curl can be performed by the standard 
spectral techniques with the use of Fast Fourier Transform. The main 
difficulty in numerical integration of the system stems from the necessity 
of transition (both direct and inverse) between the variables r and a at 
each time step. It was circumvented by the use of two independent grids 
in the r-space: a moving one (the R-grid), the motion of whose points is 
governed by (3), and a steady regular one (the r-grid), which coincides 
with the a-grid. The numerical algorithm consists of the following steps: 

(i) by integrating (3) in time, find new positions of the R-grid points; 

(ii) compute new values (1) of cj on the R-grid by finite differences; 
(Hi) by linear interpolation from the values of vorticity at nearby points 
of the R-grid, determine u> on the r-grid (for that, for each point of the 
regular grid it is necessary to find a tetrahedron, containing the point, 
whose vertices are the nearest points of the R-grid); 

(iv) solve the problem (5) to determine flow velocity v on the r-grid; 

(v) by linear interpolation, determine v on the R-grid. 
Computations are performed with the resolution of 128 3 grid points. 

In order to check numerical stability of the algorithm test runs are made 
for several initial conditions, which are ABC flows. Any ABC flow is an 
eigenfunction of the curl and hence it is a steady solution to the force-free 
Euler equation (6). They are found to remain steady in computations 
with the time step dt = 10~ 3 up to t = 4 with the relative error of 
the solution being within the 10 -14 threshold, and the Jacobian J being 
reproduced with the 10 -13 accuracy. 

An initial vorticity which we consider is a solenoidal field comprised 
of random-amplitude Fourier harmonics with an exponentially decaying 
spectrum; the decay is by 6 orders of magnitude from the first to the 
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Figure 1. The spatial minimum of (vertical axis) as a function of time (hor- 

izontal axis) at the saturated regime close to the time of collapse. Pluses show the 
values obtained in computations. 



last spherical shell in the Fourier space, the cut-off being at wavenumbcr 
8. It satisfies w/0 everywhere in the box of periodicity (this enables 
one to perform the projection (4); it is checked that this condition is not 
violated at all times during the run). This field does not possess any 
symmetry. In the course of numerical integration we monitor energy 
conservation: kinetic energy of the flow remains constant with the accu- 
racy better than 1%. For such an initial condition we observe formation 
of a peak of |u;| at a single point. At this point the Jacobian J and | ctJ | 1 
are minimal over space at all times close to the time of collapse, and the 
minimal values decrease in time to a high precision linearly (Fig. 1). (In 
this run the time step is dt = 1CP 4 for t < 0.08, and dt = 5 ■ 10~ 6 after- 
wards.) In this run the maximum of vorticity increased almost 20 times 
before integration was terminated. The final width of the peak of \co\ is 
2-3 times the length of the interval of spatial discretization (Fig. 2 shows 
a strong localization of \u>\ at the end of the run). Figure 3 illustrates 
concentration of vorticity lines and formation of a fold near the point 
of singularity. Formation of similar peaks of vorticity accompanied by 
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Figure 2. Vorticity |w| (vertical axis) as a function of the coordinates Ri and Ri 
at the plane R3 =const, containing the point of minimum of J at t = 0.08055 (close 
to the time of collapse) . 

a decrease of J to zero is also observed in several runs for other initial 
conditions of vorticity in the same class. 

At different times the global (over space) minima of J and of 
are achieved at four different points. Behavior of the Jacobian in time 
at one of these points (short-dashed line on Fig. 4) suggests that the 
second singularity can also be developing; it is not traced down to the 
time of its collapse, because this is prevented by formation of the first 
singularity. 

The peak of vorticity turns out to be narrow from the moment of 
its birth. In order to verify that it is not spurious (i.e. it emerges not 
due to a numerical instability of our algorithm) we have reproduced its 
formation in computations by a modified algorithm, with different inter- 
polation techniques employed for linear interpolation at step (iv). These 
techniques introduce some smoothing intended to inhibit formation of 
a spurious singularity. However, in the new run all numerical data has 
been reproduced with the relative precision 10~ 6 . 
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Figure 3. Isolines of restriction of the function Ri (a) on the plane oi =const through 
the point a = (77r/32; 41-7r/32; 137r/8), where the collapse occurs. Small dashes show 
downhill directons. 

To check that the given process can be considered as breaking of 
vortex lines we compute time dependencies of the Hessian of the Jacobian 
d 2 J /da a dap at the point of the minimum of J. At the final stage of the 
saturated asymptotic linear behavior of the minimum we did not find any 
essential temporal variation of its eigenvalues. This agrees qualitatively 
with the expansion (8). Figure 3 illustrates final positions of vortex 
lines near the point of collapse. Some anisotropy is observed in the 
spatial distribution of u>(r,t) near the maximum of vorticity. However, 
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Figure 4- Time (horizontal axis) dependencies of four local minima of J (a) and of 
(b). Line styles correspond to the same points on space. Solid line - the global 
minima at the final saturated regime. 

due to an apparent lack of spatial resolution we cannot claim that two 
essentially different scales emerge. The following questions also remain 
open: Why is the time of occurrence of collapse small compared to the 
turnover time? Why is the peak of vorticity quite narrow basically from 
the very moment of its appearance? 

The obtained results can be interpreted as the first evidence of the 
vortex line breaking; the collapse, which is observed numerically, is not 
related to any symmetry of the initial vorticity distribution and in par- 
ticular the collapse occurs at a single point. 
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